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d* dibaryon study is a critical test of hadron interaction models. The electro production cross 
sections of ed ed* have been calculated based on the meson exchange current model and the cross 
section around 30 degree of 1 GeV electron in the laboratory frame is about 10 nb. The implication 
of this result for the d* dibaryon search has been discussed. 
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Q\ ■ I. INTRODUCTION 

a^ : 

" ' ' The fundamental blocks of strong interaction have been confirmed to be quark and gluon. Meson, baryon, nucleus, 
O , and neutron star are the observed samples of the strong interaction matter. Theoretically it is expected that there 
2^ ' should be other samples of the strong interaction matter such as exotic quark-gluon systems, strangelet, strange star, 
and quark-gluon plasma. Experimental search for these new strong interaction matters has been undertaken for two 
' decades, even though there are some candidates of these new strong interaction matters but nothing new has been 
^SJ [ established. New facilities have been under operation or under construction for further experimental search. To help 
these search projects, new theoretical input is needed. An estimate of the d* production cross section through the 
I nd — > ird* aimed to use the LAMPF tt beam was reported in 1989 Q|. A strong production of d*, aimed to use the 
^ . existed proton machine, has been reported by Wong j^]. Intensive electron beam facilities in the GeV energy range 
are available. This paper reports an electro-production calculation of d* {IJ^ = 03+) 
, There are experimental indications of dibaryon states. A high mass 2.7 GeV) dibaryon was reported by the 
■ Sacley group Q and a low mass(2.06 GeV) one was reported by Moscow- Tuebingen group Q in addition to other 
T— I more tentative ones. H particle Q has been hunted for more than 20 years. Why is d* interesting? 
Q^ ■ Dibaryon states are closely related to the hadronic interaction, which is too complicated to be studied directly from 
the fundamental strong interaction theory, the QCD at present. Many QCD inspired models have been developed to 
. describe the hadronic interaction. 

The meson exchange model, based on meson-baryon coupling, was developed long before QCD [Q and it is still 
best at fitting the experimental data quantitatively However its validity in QCD is not clear at the moment. 
Moreover there are many phenomenological parameters fixed in the fitting process and in turn it is hard to make 
^ ] definite predictions for new physics such as the dibaryon states Q . 

Chiral perturbation effective field theory [|lO| employs Goldstone bosons, resulting from spontaneous chiral symmetry 
breaking, as the effective degree of freedom in the low energy region, and the extension to the N-N interaction is 
encouraging pl] |. Dibaryon has not been studied in this approach yet. 
$^ ' L.Ya. Glozman, D.O. Riska and G.E. Brown proposed that the Goldstone boson is not only an effective degree 
of freedom for describing hadronic interactions but also good for analyzing baryon internal structure. In their model, 
constituent quarks and Goldstone bosons are used as the effective degrees of freedom. Up to now the application is 
mainly restricted in the baryon spectroscopy except a study on the origin of the N-N repulsive core. 

A. Manohar and H. Georgi [^, however, have argued that between the chiral symmetry breaking scale(~ 1 GeV) 
and the confinement scale(~ 0.2 GeV), the effective degrees of freedom are Goldstone bosons, constituent quarks and 
gluons. Such a hybrid quark-gluon- meson exchange model has been developed to describe nuclcon-baryon interactions 
and a semi-quantitative fit has been obtained . Some dibaryon states have been studied with this model El] 
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A constituent quark and effective one gluon exchange model Q describes hadron spectroscopy quite well |17|, but 
only the repulsive core of the N-N interaction was obtained when this model was extended to study hadron interactions 

®- . . rr. 

The MIT bag model uses the current quark and gluon to describe hadron internal structure |19| |. It was used 

extensively in the study of dibaryon states in the early 1980's and resulted in an explosion of dibaryon states [ pO| . 
It was realized latter that the unphysical boundary condition should be modified IM. One modified version is the 
R-matrix method and the other one is the compound quark model approach |]23t] . 

The Skyrme model ||2^ has also been used to study hadron interactions |^ and dibaryons |^6| . Additional models 
might exist that should be added to the list. 
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It seems hard to discriminate these models just by the hadron spectroscopy and the existed scattering data of 
hadron interactions. TheoreticaUy it is also hard to justify which effective degrees of freedom are the proper ones. 
On the other hand the well known phenomena, that the nucleus is a collection of nucleons rather than quarks, and 
that the nuclear force has similarities to the molecular force except for energy and length scale differences, have not 
been explained by any of these model approaches. 

A pure quark-gluon model description of the N-N interaction has been developed . It starts from a multiquark 
system and demonstrates that in the N-N channels, it is energetically favorable for the system to cluster into two 
nucleons and that the nuclear intermediate range attraction is caused by quark delocalization similar to the electron 
delocalization which induces intermediate range molecular attraction. In the 9 and 12 quark systems with quantum 
numbers of ^H, ^He and '^He the nucleon clustering has been verified as well as the two nucleon system [g8[. This 
model has been extended to N-k and iV-E interactions p9| and the results show that the quarks delocalize properly 
in different channels to induce qualitatively correct N-NXh = 10, 01, 11, 00), A^-A (JI = 1^, 0^), and N-T, (JI = 
li, Oi, If, 0|) interactions. For other channels §], such as AA (JI = 30), Arj(3|), AE*(3i) and AS*(31), it is 
energetically favorable for the quarks to merge into quark matter instead of two baryons, and there are often strong 
effective attractions in these channels. In the A A (JI = 30), the d* channel, the effective attraction is so strong that 
the total energy of the system(~ 2.1 GeV) is near the NNtt threshold; therefore the d*(JI = 30) might be a narrow 
resonance state [||. 

Different model approaches give quite different mass of d* . The meson baryon coupling model 1^ gave a A — A 
binding of 11 — 340 MeV depending on the coupling constants of p and uj and the hard core radii(0.2 — 0.3 fm). If the 
hard core radii are larger than 0.38 — 0.40 /m,the binding energy will be less than 10 MeV . The quark-meson-gluon 
hybrid model obtained a binding between 37 MeV and 80 MeV The pure gluon exchange model got 18 MeV 
binding even though there is a repulsive core in many other baryon-baryon channels |31|. MIT bag model will give a 
d* mass of 2340 MeV if the bag radius is adjusted to give the minimum. On the other hand the R-matrix version of 
the modified bag model give a mass as high as 2840 MeV [Q. The skyrmion model obtained a very weak binding 
~ 10 MeV 1^^. Therefore d* study will provide a critical test of hadron interaction models. (This was already 
pointed out in the US Long Range Plan for Nuclear Science in 1996 [3^].) 

Quark delocalization plays a vital role in lowering the d* mass. In a variational calculation this is nothing else 
but just a method to enlarge the variational Hilbert space. Therefore it might be a general property of quantum 
mechanics. If it is really realized in the nucleus (nucleon swollen explaination of the EMC effect |3J] might be taken 
as an evidence of quark delocalization), it has far-reaching implications however. The phase transition from nuclear 
matter to quark-gluon plasma would be at best a second order one or just a crossover as happened in the transition 
from atomic gases into plasma and would make the already hard identification of this phase transition even harder 
B . d* search is a critical test of the quark delocalization mechanism. 



II. ELECTO-PRODUCTION MECHANISM OF D* 



d* is a spin 3 state. Its dominant hadronic component is AA. To produce a d* from a nuclear target, one has to 
change two nucleons into two A's. A virtual photon exchange can only excite one nucleon into a A. Double photon 
exchange is a fourth order QED process. Its contribution to the cross section will be proportional to a'', a too small 
effect. Here a = e^/ {he). Therefore the electro- production of d* is critically dependent on the gluon exchange current 
in the quark-gluon description and on the meson exchange current in the hadronic description. The gluon exchange 
current is unique as depicted in Fig.l. 
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FIG. 1. Quark gluon exchange current 




FIG. 2. Quark meson exchange current 

However as mentioned before that some models emphasized the Goldstone boson quark couphng. Then there will 
be meson exchange current as depicted in Fig. 2 even within a baryon. This makes the exchange current in the quark- 
gluon description rather model dependent. J.A.Gomez Tejedor and E. Oset(GO) calculated the ed e'pA and 70? 
AA prni^TT" cross sections with the hadronic degree of freedom ||3^. This approach is complicated by the many 
effective meson-baryon couplings. However those coupling constants have been fixed by the experimental data. We 
will take GO approach to do the d* electro-production cross section calculation. Based on the isospin and angular 
momentum conservation and taking into account the results of GO, the following Feynman diagrams are included in 
our calculation, 
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A- A++ 



FIG. 3. KroU-Ruderman process 
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FIG. 4. Meson exchange current 
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FIG. 5. NN intermediate state 
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FIG. 6. NN* intermediate state 



The d* is an isospin 1=0 state, 

I d*) = i(A++A- - A+A° + A°A+ - A-A++). (1) 

The KroU-Rudcrman term Fig. 3 does not contribute due to the cancehation of the 7r+ and tt^ exchange between 
different components of d* and d. The meson exchange term Fig. 4. does not contribute for the same reason. Therefore 
only Fig. 5 and 6 contribute to the d* production. 

For the NN* intermediate state, only N*(1520 /J^ = M^) been included. Because the results of GO |^ show 
that the contribution of NN*(1520) might be more important, the NN*(1440) term will be left for further refinement. 
The AA and the D-wave components of deuteron are neglected temporary. 



III. MESON EXCHANGE CURRENT AND CROSS SECTION 



The general electron scattering cross section formula of Donnelly and Raskin will be used to calculate the 
inelastic d* production. 



^ = aMi>^LR% + ^tRI + ^TtR"}! + VTLR)t)frel (2) 



where ctm is the Mott scattering cross section, 
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2\ 2 

q 



2 \q- 



l^TT = -77 (^), (6) 



1 (t-\J-(t 



R% = \piq)M^ (8) 

Rl^\J{q,+l)M^ + \Jiq,-l)f^f, (9) 
= 2Re{r{q, +1)/. J(g, -1)/.}, (10) 
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= -2Re{p*iq)f,iJiq, +1);, - J(g, -1)/.)}, (11) 

/.ec = l + ^sin^|. (12) 

free is the recoil correction. The four vector current J^{q)fi is the Fourier transformed four vector transition current 
of the target, 



J''{q)H = j d're^^-%f\J^{r)\i), (13) 
J{<l)r^^ Ji<l,m)eiq;l,m), (14) 



m=0,±l 

1 



e((f;l,0) = U2, e((f; 1,±1) = =F — ("x ± iwy), (15) 



q X (k X k') ^ k X k' ^ q , „, 

Ux = — , Uy = — —, Uz = -r-r, 16 

\qx{kxk')\ " \kxk'\ \q\ 
p{q)f,= ^^J{q,Q)f,. (17) 

The last Eq. is obtained from the four vector current conservation, g, fc, and k' are the four momentum transfer, the 
initial and final momentum of the scattered electron depicted in Fig. 7. q = k — k' . 9e is the electron scattering angle 
in the lab system. Mt is the mass of the target, mass of deuteron in our case, e (e') is the initial(final) energy of the 
scattered electron, and w = e — e'. 




FIG. 7. The electron-deuteron inelastic scattering 



For the unpolarized electron scattering, the cross terms Rj^ and Rj^" do not contribute. The angular momentum 
and parity conservation {IJp = 01+ for d and 03+ for (i*)restrict the multipole moments further. Only the following 
5 multipole moments contribute to the ed — s- ed* process: 



fFM) 

Ffim 



/47r , , 
iAtt 



|M2(k1)||l>, 

|M4(k1)|ll), 

|Tf(k1)l|i>- 



(18) 
(19) 
(20) 
(21) 
(22) 
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The multipole moments are defined through the following multipole expansion, 

Amn = ^v^$^\/2JTTi''(/|M^o(|9l)|i), 

J>0 

J{q,±l)si = -s/2^Y.^^^i' ((/m±i(l3l)N) ± {f\Ti:±Ami)) , 

J>0 

MjM) = J dVMj„(|glr)p(r), 

T^rM)^ Jd'rMji^m{W)-Jir), 
Mjm{W)=jjmr)Yj^{f), 

MjLmiW) ^ Jj{\^r)YjLm{f), 

YjLr^{r) = [YL{f)^e{q;l)]j^. 
The transition current resulted from Fig. 5(a) is 

({q + q'f) J , rt . Tl6^ {q+ p, - p,. ) , 

A similar transition current resulted from Fig. 6 (a) is 

((9 + Q') ) 7 ^1 . T^6' {q+p, - p,, ) 

" ~ 1 / /* 



Ul + qf) -—77^ . T^6' {q+ p, - p,, ) . 

^ [q + q ) - fj- + «e 



Fig. 5(b) and Fig. 6(b) will give similar transition currents. 

In obtaining these transition currents,the following efi'ective Lagrangian have been used, 

LNAn = -—■^isj (a^c/.^) r^^*iv + h.c. 

i^N'fN'A^4'^T^''^A + h.C. 



In these expressions <J>, ^n, ^a, ^n*, and stand for the pion, nucleon, A, A''*(1520), and photon fields 
respectively; Mjy and ^ are the nucleon and pion masses; a and r are the usual i spin and isospin Pauli operators; 

and are the transition spin and isospin operators from ^ to | with the normalization, 

d,M\Sl\l,m) =C (l,l,^;m,:^,M) (39) 



2' ' V2 2 

All the effective coupling constants, form factors, and the propagators are taken from GO |35j. They are copied 
here to make this report self contained and can be read without further check of many references. 
1.79 proton 
-1.91 neutron 



/* - 2.13 /at. A. = 0.677 

. _ f 0.108 proton 1 - - / -0-049 proton 

^ I -0.129 neutron | 9'^ - \ _o.0073 neutron 



For the nucleon propagator, only the positive energy part is retained where E{p) = \J Mf^ + ; for the N* , the 
finite width T has been included where 

5 = _ Y{VS) = r(M^)gL(V^)/'zL(Miv*). (40) 

The form factors are taken from GO |35j to keep the model consistent. 

^^(<?') = -ri h K^l2bQMeV (41) 



Sachs's form factors are given by 



Gf(,^)^^^-^ (42) 



with A2 = 0.71 GeV^; /Zp = 2.793; /i„ = -1.913. 

The relation between F^{q'^) (Dirac's form factor) and G^^{q^) is : 



1 tJ-p ' 



F?iq')^GUq')^-^—^ (43) 



and Fl' = 0. 

To calculate the hadronic d d* transition current, a single Gaussian wave function with a size parameter 
b* = 0.7 fm is used to approximate the d* internal motion. For the deuteron d a three Gaussian fitted to the 
ground state properties has been used, three size parameters are b = (4.0975, 1.8252, 0.8837) fm with normalization 
coefficient c = (0.31491, 0.49716, 0.36926). The d* mass is taken to be 2.1 GeV. Due to the special spin(l for d 
and 3 for d*) and orbital angular momentum(both 0) internal structure and the spin property of the meson exchange 
current, only F2 , F^ and F^ three transition form factors remain after the integration over the internal spin and 
orbital variables of d and d* . Fig. 5 contributes one F2 , one F^'^ and two _F^(from convective and magnetic current 
respectively). Fig. 6 contributes one and one F^ term. 



IV. RESULTS AND DISCUSSIONS 



The calculated d* electro production cross sections are shown in Fig. 8 . The four curves correspond to four electron 
energies 0.8, 1.0, 1.2 and 1.5 GeV. The shape of the differential cross section is dominated by the Mott cross section 
but modulated by the inelastic transition form factors. The value around 30 degree in the laboratory frame is about 10 
nb for 1 GeV electron. In Fig. 9, 10 and 11, the transition form factors Ff, F2 and F*^ are shown. The dashed curves 
correspond to the contribution of NN intermediate state while the full curves correspond to the sum of contributions 
of two intermediate states. The NN intermediate state is the dominant one except at large angles, where the NN* 
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contribution is more important. Fig. 12 shows the differential cross section due to the NN intermediate state only as 
well as a comparison to the full one. 

The produced d* will decay into NN and NNtt. The cross section, 10 nb around 30 degree for 1 GeV electron, is 
about two order smaller than that of quasi elastic ed e'NN and A resonance production ed e'NA e'NNir 
processes. Therefore the normal measurement of the inelastic electron scattering can not find the signal of d* even it 
is existed, i.e., the d* signal will be buried in the quasi elastic or A resonance background. A special kinematics and 
detector system must be studied further both theoretically and experimentally in order to pin down the weak signal 
of d* resonance from the strong background. 




e(deg) 



FIG. 8. ed — » ed* differential cross sections 
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FIG. 12. Same as FIG. 9 for differential cross sections da/dQ 



In our model calculation, some simplifications have been assumed. The D wave component of deuteron has been 
neglected. Even though it is a small component, its contribution to the d* production might not be small enough to 
be neglected, because only a recoupling of spin and orbital angular momentum will be able to transit the deuteron 
"^Di into a ^D^ NN component of d*. The initial state correlation, i.e., the A — A component of deuteron is an 
even smaller component, but only a spin recoupling is needed to transit it into d*, therefore its contribution should 
be checked as well, d* is a six quark state y] rather than a A — A bound state. To model it as a pure two A 
bound state and use a single Gaussian wave function to describe its internal structure will overestimate the calculated 
cross section. The contribution of A^A^*(1440) intermediate state should be added especially in the large angle part. 
Due to these approximations the calculated cross section is an estimate of the d* electro production. On the other 
hand all of the above mentioned corrections are minor effects, the order(10 nb around 30 degree for 1 GeV electron) 
obtained from this calculation might be stable against these fine tunes. Further calculation is going on, especially a 
quark model calculation is also doing to check if a single Gaussian approximation of the d* internal wave function has 
overestimated the cross section.. The results will be reported later. 

The electro production of d' dibaryon is being measured. This calculation can be extended to that process and will 
be a good check not only of the electro production model used here but also of the d' dibaryon analysis itself. 

Very helpful discussions with C.W .Wong, T. Goldman and Stan Yen are acknowledged. We are also greatly 
indebted to J. A. Gomez Tejedor and E. Oset for their helpful private communications. 

This research is supported by NSF, SSTD and the post Dr. foundation of SED of China. Part of the numerical 
calculation is done on the SGI Origin 2000 in the lab of computational condensed matter physics. 
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